Chiral corrections to the scalar form factor in B n — ► D„ transitions 
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We consider chiral loop corrections to the scalar form factor in B q —> D q lvi decays. First we 
consider chiral corrections to the 1/tuq suppressed operators and then we propose the procedure for 
the extraction of the relevant form factor using lattice QCD results. In the case of B s — » D a lui decay 
we find that effects of kinematics and chiral corrections tend to cancel for the scalar form factor 
contributions. In particular the 1/mg suppression of chiral corrections is compensated by the 517(3) 
flavor symmetry breaking corrections which can be as large as 30%. The calculated corrections are 
relevant for the precise determination of possible new physics effects in B q — > D q lvi decays. 

PACS numbers: 12.39.Fe, 12.39.Hg, 13.20.He 
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I. INTRODUCTION 

The b — > c transitions have been traditionally impor- 
tant for the determination of the V c b CKM matrix ele- 
ment. For the exclusive determination at the B facto- 
ries, the B — > D* ev decay channel is the most suitable. 
In addition to being the dominating B meson scmilcp- 
tonic decay mode, heavy quark symmetry (Luke's theo- 
rem) protects the normalization of the only relevant form 
factor at zero recoil point from he leading 1/ttiq power 
corrections [1(. Furthermore, its shape can be extracted 
from the experimentally measured decay spectrum. How- 
ever, with the modulus of V c b determined to a percent 
level, the focus is shifting toward constraining possible 
new physics contributions also in charged current transi- 
tions - traditionally a domain of neutral current studies. 
B — > Dtv is one such promising process, where the large 
tau mass makes it sensitive to hclicity suppressed con- 
tributions as such coming from the exchange of charged 
scalars @, S 0, Q • This calls for high precision theoret- 
ical estimates for both the dominating vector (G(w)) as 
well as the relative scalar form factor contribution (A(w)) 
describing the non-perturbative QCD dynamics of this 
transition. First such estimates on the lattice have al- 
ready been performed p, B S Q ■ 

At the same time, high luminosity hadron machines 
(Tevatron and LHC) will take the central stage in the 
coming years also in the area of precision heavy flavor 
studies. Their higher energies allow to study processes, 
unavailable at the B factories. B s — * D s transitions 
are examples of these. Although to our knowledge no 
non-perturbative lattice estimates exist at the moment 
for the required form factors, they are actually better 
suited for lattice studies than B — > D transitions, as 
large chiral extrapolations to the physical light valence 
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quark masses and the associated systematical errors can 
be largely avoided. 

In this note we estimate the leading chiral symmetry 
breaking corrections, governing the differences between 
the vector and scalar form factors in B s — > D s and 
B — + D transitions. These corrections can on the one 
hand be used to guide lattice QCD studies in their chiral 
extrapolations. On the other hand we use them to esti- 
mate qualitatively the relative scalar form factor values 
in B s — ► D s transitions. 

The leading order (LO) SU(2>) chiral corrections at 
leading and next to leading order (NLO) in 1/ttiq to 

the B q — > D q ' scmileptonic form factors have previously 
been computed in [lj|. However phenomenological dis- 
cussion at the time was limited by the lack of experimen- 
tal and lattice information available on the sub-leading 
form factor contributions. We are now able to provide 
the first more reliable estimates of these contributions 
and also study their impact on B s — ► D s tv phenomenol- 
ogy 



II. FRAMEWORK 

A. 1/mQ power expansion of heavy quark currents 

The matrix element of the vector b — ► c quark current 
between a B and a D meson of velocity v and v' respec- 
tively, can be parametrized in terms of two form factors 

(D(v')\&y^,b\B(v)) = y/m B m D [h+(w)(v +v% 

+h-(w)(v-v%] , 

(1) 

where w = v ■ v'. In heavy quark expansion they take the 
form [n[ 

h+(w) = £H+0(e Q ), (2a) 

h-(w) - (e c -e b )[-M(w)+2&(w)]+0(e 2 Q ), 

(2b) 



where eq = l/2m,Q and A = m# Q — rriQ. The £,(w) 
and £3(11)) form factors are defined through the match- 
ing of the LO and NLO quark currents in HQET to 
HMxPT [HEH. For the leading (£(iu)) form factor (the 
Isgur-Wise function) one gets 



h'Vh 



-£,{w)Tr[H'TH] 



(3) 



where h(x) = cxp (—irnqv ■ x)P + Q(x) is the HQET 
heavy quark field, P+ = ('j-v+l)/2 is the HQET velocity- 
spin projection operator, while H(x) = P + [P*(x) ■ 7 — 
P{x)^ 5 ] is the HMxPT l/2~ field of vector and pscu- 
doscalar mesons. All the primed (') quantities refer to 
the primed velocity v' and flavor (mass itiqi). The lead- 
ing order contribution to the h-(w) form factor origi- 
nates from the matching of QCD to NLO HQET cur- 
rents. At tree level, the Lorentz structures of T are the 
same on both sides and we are interested in the case 

r = K_ A = 7^(1 - 75) M 

(H'\hTj a W a h\H) = Tr[HTP-j a H£ a (v,v')], 

(4a) 
(H'\(-W a )h'j a Th\H) = Tr[t a (v',v)H'~f a P'_TH}, 

(4b) 

where 

£ a (v,v') = £+{w)(v + v') a + £-(v - v') a - £ 3 (w)~f a , (5) 

and £ a = 7o£i7o- The following identities can be proven: 
e_H = H(w)/2, t + (w)(w + l)-(w-l)Z-(w)+Z 3 (w) = 
0, reducing the number of independent form factors to 
the two (£(w), £,3(w)) in eq. J2b|. Additional 1/rriQ form 
factors originate from the insertions of NLO HQET in- 
teraction Lagrangian. They contribute however only to 
h+(w) at 0(eq) and as such do not interfere in the LO 
determination of h-(w). 



B. Chiral expansion in HM^PT 

The framework of the chiral heavy meson Lagrangians 
and procedures involved in the computation of loop cor- 
rections in the chiral expansion is by now standard and 
has been explained previously (c.f. [3|)- Novel however 
is the consideration of the HMxPT Feynman rules for 
the sub-leading current operators from the previous sec- 
tion [22|]. Since all of them involve power corrections 
to both initial and final states, their linear combinations 
entering expressions (|2al2b[) in general mix under chiral 
corrections. The diagrams needed to be computed are 
the self-energy contributions for both the initial and fi- 
nal states as well as the "sunrise" diagrams with pseudo- 
Goldstone exchanges between the initial and final states 
connected to the weak current operators (fig. [I]) . At LO 
in the chiral expansion no additional pscudo-Goldstonc 
fields can be emitted from the weak current operator 
vertices. In the computation of the power suppressed 
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FIG. 1: Weak vertex correction "sunrise" diagram. Crossed 
box represents effective weak current operator vertex, while 
filled black circles represent effective strong vertices. 



chiral corrections to h + {w), one would need to con- 
sider contributions from the NLO HMxPT interaction 
Lagrangian [lfj. These are however further suppressed 
when considering LO chiral corrections to the sub-leading 
form factor h-{w). 

In our phenomenological discussion we consider two 
cases: if one is to use the one-loop formulae to extrapo- 
late lattice QCD results on the scalar form factor to the 
chiral limit, one has to use the SU(2) theory. Namely, 
it has been shown previously that LO chiral symmetry 
breaking contributions in an SU(3) limit cannot be reli- 
ably disentangled from the contributions of the lowest ly- 
ing 1/2+ excited heavy meson states [1J, [jji [l6[ . On the 
other hand with a number of additional assumptions one 
can provide a qualitative estimate of the B s — » D s sub- 
leading form factors even before lattice QCD results on 
these become available using SU(3) theory predictions. 



III. PHENOMENOLOGY 



The B, 
in terms of w — vb ■ Vd as 

dT(B -► D£v) 



q —r D q lv partial rate in the SM can be written 

G 2 F \V cb \ 2 m% 

a pv{w) (6) 



dw 
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where t(w) = m B + m 2 D — 2wm^mB and we have decom- 
posed the rate into the vector and scalar Dalitz density 
contributions 
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Possible new charged current scalar interactions only 
contribute to pg, introducing a shift ps(w) — ► [1 + 
t(w)5pfp]ps(w), where Spjp parametrizes the new physics 
contribution @, [j, Q- The G(w) and A(w) form factors 



can be written in terms of h+(w) and h-(w) as 



G(w) 

A(w) 



h+(w) — 



m B - m D 



h-{w) 



rn D + tub 
iriB — ttld w — \ h_(w) 



n%B + Tod w + 1 h+(w) 
tub — triD h—(w) 
tub + mD h+(w) 



(9) 



(10) 



The benefit of such form factor decomposition is also that 
G(w) completely dominates B — ► Dtv decay rates with 
t = e, n and its shape can actually be extracted from the 
corresponding decay spectra. On the other hand, A(u>) 
is only relevant for B — > Dry and in addition depends 
on the non-perturbative QCD dynamics only through the 
ratio of form factors h-/h+, making it highly suitable for 
lattice QCD studies. 

G(w) receives chiral corrections already at LO in 1/toq 
(the well known chiral corrections to the Isgur-Wisc func- 
tion £(w) [171 ll8J) making the sub-leading contributions 
perhaps less important. Only at zero recoil (w = 1), 
Luke's theorem protects h+(w) from 1/toq corrections 
(£(1) = 1 by construction and thus invariant under chi- 
ral corrections), and here the leading order chiral cor- 
rections to G(l) actually come from 1/toq suppressed 
h-(w) contributions. The case is similar for A (it;), whose 
leading order value is strictly determined by kinematics 
throughout the kinematic region. Here the chiral correc- 
tions to sub-leading terms (proportional to h_/h + ) are 
dominating the chiral extrapolation of lattice studies as 
well. Both cases underline the importance of controlling 
the chiral corrections of the h-(w) form factor. We turn 
to this issue first. 



current operators, needed to cancel the UV divergences of 
chiral loops. We have used eqs. (|2al2b[) to rewrite £(3) (w) 
in terms of ti± ee (w), which is valid up to 1/toq correc- 
tions. Since r(l) = l. chiral corrections to /i+(l) at LO 
(and also at NLO [lj|) in 1/toq vanish as they should, 
and the dominating contributions come at 1/toq [191 ] - 
The value of h-(l) however is not protected by heavy 
quark symmetry. Note also that the one-loop chiral cor- 
rections to h- (w) are proportional to a different combi- 
nation of £3(11)) and £(w) than the tree level value result- 
ing in an admixture of h + (w) in the /z of m^ running of 
h—(w). This can be most easily understood at zero recoil, 
where the value of £(1) admixture in ft-(l) is protected 
from receiving chiral corrections by heavy quark symme- 
try. Therefore they should be proportional exclusively to 
£3(1), which is the case. Actually the eigenvectors of the 
two coupled equations (|llalllb[l in the whole kinematic 
region are precisely £,(w) and £3 (w). For completeness, 
we give their explicit chiral LO corrections formulae in 
the appendix. Incidentally, when evaluating the A(w) 
form factor, only the ratio h~/h+ is needed, whose chi- 
ral corrections take on a simpler form, as can be easily 
inferred by dividing eqs. (|1 la,|l lb[) and then expanding 
up to linear terms in m„ log m\ 
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A. Lattice QCD extrapolation formulae 

Below we give the leading logarithmic chiral correc- 
tions relevant for lattice QCD extraction of the h+(w) 
and h-(w) form factors in the Nf = 2 theory (r(x) = 
log(a; + Vx 2 - l)/Vx 2 - 1) 
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hi* 



-Tree 



l + 3.g 



1 - 3g 



2 r{w)-l 2 



(4^/.) 2 



2 + Y+[r(w) + l] 



where Y + = (e c — £f ) )A/i^ ree //iZ ¥ee and the w dependence is 
implicit. Here and in the rest of the text Ci(^(, w) denote 
the sums of local analytic counter-terms, which cancel 
the fi dependence of the chiral log pieces. As stressed 
by the notation, they will in general have a non-trivial 
u> dependence, originating both from the finite analytic 
residuals of chiral loops as well as from local NLO chiral 



Again one can understand the form of chiral correc- 
tions by rewriting the ratio h—/h+ — > (e c — e&)(— A + 
2&( w ) / £,( w )) ■ Since A is defined independently of 
spectator effects, it should not receive chiral correc- 
tions. Therefore they should again be proportional to 
£3(11}) /£(w) throughout the kinematic region, as can be 
readily checked by using eqs. (|2al2bj) . 

The same sub-leading IW functions enter also other 
b — > c hadronic matrix elements. In particular, the A(e c — 
£b) contribution introducing explicit heavy quark mass 
dependence or equivalently, the mixing of leading and 
sublcading form factor contributions can be avoided by 
rewriting the chiral extrapolation formula (jllbp in terms 
of the tree level matrix element (D*(e\v')\T\B*(e,v)) 
part proportional to — e ■ e'(v — v') - we denote it with 
h*_ [23]. In heavy quark expansion it can be written as 
h*_(w) = -A(cc - eb)t(w) + C(e Q ) [ij. Immediately 
2 c_(tt, w) \ , (lib) w e see, that its leading chiral corrections should coincide 
with those of h+ as can also be checked by computing the 
two sunrise diagrams contributing to this vertex correc- 
tion. The chiral extrapolation equation for h- (w) (jllbp 
on the other hand remains of the same form with the 
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replacement -Y + -> Y* = {h*_/h_) 1Ycc . The goal of lat- 
tice simulations is then to determine h^ ee (w) or alterna- 
tively h_ ee (w), and Ci(w) by fitting the above formulae 
to the numerical values of h±({m 7r })(w), /il({m 7r })(w) 
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FIG. 2: Estimate of £3(10) from lattice data on h±(w) [a, Sl- 



at several values of m, for each considered to. Using 
eqs. (9) and (10) (and (7) and (8)) one can thus obtain 
all physically relevant hadronic variables needed in the 
search for possible charged current scalar interactions in 

B q — > D q TU. 



B. B s — > D s estimate for A(w) 

We can revert the game and use the already available 
information on the h-(w) form factor in B — > D, to esti- 
mate the size of the sub- leading form factor £3 (to) . Then 
we can use this information together with a number of 
assumptions to estimate A(to) in B s — ► D s . Namely, 
using eqs. (|2al2b|) with A ~ 0.5 GeV [H] and the MS 
charm and b quark masses m c (m c ) = 1.2 GeV, mb(rrib) = 
4.2 GeV [21j we can extract from available lattice data on 
h + (w) and h-(w) [H, Q the size of £(10) at leading order 
in 1/rriQ and consequently £3(11)) for the calculated lattice 
points as shown in figure O Then we use the computed 
chiral corrections to £(to) and £3(10) in the Nf = 3 limit, 
neglect the analytic terms at a fixed scale of /i = 1 GeV 
and also neglect contributions from excited heavy meson 
states to obtain an estimate for A(w) in B s — > D s . The 
dominating contribution naturally comes from the LO 
kinematic factor since mr> s jmB s 7^ Tod/ tub- This con- 
tribution amounts to roughly 3% supression from unity in 
the ratio A Bs ^ Ds (to) / A B ^ D (w). Additional chiral con- 
tributions come from the chiral corrections to the HQET 
form factors entering A (to) and can be put into the ap- 
proximate form 

h'-h* = 1 + 2g2 (l-n)(l + rH) 



hih 8 , 
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x -rrvlog- 
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1 2, m r, 

2 m " l0g ^ 



C.t. , 

(13) 



where I/I5 = 1 — 2^3/A^ and c.t. denotes additional 
counter-term contributions. Both kinematic and chiral 
log effects arc shown in fig. [3] [24j . Only uncertainties 



FIG. 3: Estimate of A Bs ^ Da (w)/A B ^ D (w) using Nf = 3 
logarithmic chiral corrections and information from lattice 
data on G(w) and A(w). The leading order kinematic ef- 
fect is shown in black, while chirally corrected result is shown 
in red. Only uncertainties coming from lattice inputs [a @| 
are taken into account in the error bars displayed. 



coming from lattice inputs are taken into account in the 
error bars displayed. We see that the effects arc comfort- 
ably small, point at suppression smaller than 5%, and 
are consistent with unity ratio within the current large 
errors throughtout the available to region. The chiral cor- 
rections are supressed by (2£ 3 — A£)(e c — e&) ~ 0.1. How- 
ever, this is compensated by the large SU(3) breaking 
corrections which are of the order 30%. A more reliable 
estimate is of course expected from the lattice. 



IV. CONCLUSIONS 

Motivated by the possible sensitivity to the presence 
of new physics in B q — > D q rv decays 0, 0, 0, Q we 
have investigated the effects of chiral corrections to the 
relevant (scalar) form factor. The explicit chiral behavior 
of the computed chiral corrections can be used to guide 
future lattice computations in approaching the physical 
regime for the light quark masses. 

Based on results of our previous studies [lj, [l5|, [lj| we 
calculate the leading chiral logarithmic behavior of the 
h-(w) and h + (w) form factors, by including only pionic 
contributions. We consider chiral corrections of the op- 
erators appearing at 1/ttiq in matching QCD to HQET 
and which contribute to the scalar form factor. We pro- 
pose to determine the matrix elements for the B — > D 
and B* —* D* transitions with the same kinematics de- 
pendence (v — v') on the lattice simultaneously. In this 
way one can extract the scalar form factor at tree level in 
the chiral expansion, as well as the relevant counter-term 
needed to make predictions in the physical regime. 

Finally we consider chiral corrections to the scalar form 
factor in the B s — ► D s transition and conclude that 
SU(3) flavor symmetry breaking effects are significant 
(up to 30%) and can thus partly compensate the 1/uiq 
suppression of chiral corrections. They are comparable 
in size with kinematical effects but tend to be of the 



opposite sign. This points towards longitudinal/scalar 
contributions very similar to B — ► D transitions. 



matrices and a summation of the repeated index i is as- 
sumed. As in the text r(x) = log(a; + yjx 2 — l)/\/x 2 — 1. 
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APPENDIX A: SU(N) CHIRAL CORRECTIONS 
FORMULAE FOR £(w) AND £ 3 (w) 

Below we give the complete expressions of the chiral 
1-loop corrected form factors £(w) and £,3(11)), which are 
eigenvectors of the chiral running and are needed for the 
LO lattice chiral extrapolation of h±(w). They are given 
in a SU(N) theory where U denotes the SU(N) generator 
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Here c.t. denotes the corresponding counter-terms, which 
cancel the jj, dependence of the log terms. The counter- 
terms are analytic in rrii and will generally themselves 
have a non-trivial w dependence, so they should be fitted 
for each value of w independently. This is also true for 
£3(11)) despite the fact, that its LO non-analytic chiral 
corrections are w independent. 
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